Astronomy & Astrophysics manuscript no. triplc2 


February 2, 2008 


(DOI: will be inserted by hand later) 





The three-point correlation function of cosmic shear. II: Relation 
to the bispectrum of the projected mass density and generalized 

third-order aperture measures 

Peter Schneider 1 , Martin Kilbinger 1 , and Marco Lombardi 1 ' 2 

1 Institut f. Astrophysik u. Extr. Forschung, Universitat Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 
e-mail: peter@astro.uni-bonn.de 

2 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85741 Garching, Germany 

Received ; accepted 

Abstract. Cosmic shear, the distortion of images of high-redshift sources by the intervening inhomogeneous 
matter distribution in the Universe, has become one of the essential tools for observational cosmology since it was 
first measured in 2000. Since then, several surveys have been conducted and analyzed in terms of second-order 
shear statistics. Current surveys are on the verge of providing useful measurements of third-order shear statistics, 
and ongoing and future surveys will provide accurate measurements of the shear three-point correlation function 
which contains essential information about the non-Gaussian properties of the cosmic matter distribution. 

We study the relation of the three-point cosmic shear statistics to the third-order statistical properties of the 
underlying convergence, expressed in terms of its bispectrum. Explicit relations for the natural components of 
the shear three-point correlation function (which we defined in an earlier paper) in terms of the bispectrum are 
derived. The behavior of the correlation function under parity transformation is obtained and found to agree with 
previous results. We find that in contrast to the two-point shear correlation function, the three-point function 
at a given angular scale 6 is not affected by power in the bispectrum on much larger scales. These relations are 
then inverted to obtain the bispectrum in terms of the three-point shear correlator; two different expressions, 
corresponding to different natural components of the shear correlator, are obtained and can be used to separate E 
and B-mode shear contributions. These relations allow us to explicitly show that correlations containing an odd 
power of B-mode shear vanish for parity-symmetric fields. Generalizing a recent result by Jarvis et al., we derive 
expressions for the third-order aperture measures, employing multiple angular scales, in terms of the (natural 
components of the) three-point shear correlator and show that they contain essentially all the information about 
the underlying bispectrum. We discuss the many useful features these (generalized) aperture measures have that 
make them convenient for future analyses of the skewness of the cosmic shear field (and any other polar field, 
such as the polarization of the Cosmic Microwave Background). 
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1. Introduction 

Recent surveys have measured second-order cosmic shear statistics with high accuracy, owing to the large sky area 
covered, and thus the large number of faint galaxy images (e.g., van Waerbeke et al. 2001, 2002; Jarvis et al. 2003a; 
Hoekstra et al. 2002). With surveys of this size, it now becomes feasible to obtain higher-order cosmic shear statistics 
which probe the non-Gaussian features of the cosmic shear field. These higher-order statistics are particularly useful 
in breaking near-degeneracies in cosmological parameters which are present at the level of second-order statistics. 
Bernardeau et al. (1997), van Waerbeke et al. (1999) and others pointed out that the skewness of the convergence 
underlying the cosmic shear field can break the degeneracy between the density parameter f2 m and the normalisation of 
the matter power spectrum, expressed in terms of the r.m.s. density fluctuations as on a scale of 8/i -1 Mpc. However, 
the convergence cannot be observed directly, but needs to be inferred from the observed galaxy image ellipticities 
which yield an estimate of the local shear. The dispersion of the shear in a (circular) aperture, frequently used to 
quantify second-order shear statistics, cannot be generalized to a third-order statistics. Schneider et al. (1998; hereafter 
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SvWJK) have defined an alternative cosmic shear measure, the aperture mass, which is a scalar quantity that can be 
directly obtained from the shear, and therefore is particularly suited to define higher-order statistics. 

Recently, interest in higher-order cosmic shear statistics has been revived. The three-point correlation function 
(3PCF henceforth) of the shear contains all the information on the third-order statistical properties of the shear field, 
and therefore is of prime interest. In addition, it can be obtained directly from the observed image ellipticities and, in 
contrast to the aperture mass statistics, is insensitive to holes and gaps in the data field. However, the shear 3PCF is 
a function with 2x2x2 = 8 components (since each shear has two independent components) and 3 variables (e.g., 
the sides of the triangle formed by the three points) and therefore difficult to handle. Bernardeau et al. (2003) defined 
a specific integral over the 3PCF and applied that to the VIRMOS-DESCART survey in Bernardeau et al. (2002) to 
obtain the first detection of a non-zero third-order cosmic shear signal. Using the same observational data, Pen et al. 
(2003) calculated the skewness of the aperture mass, where the latter has been obtained from integrating the shear 
3PCF. Jarvis et al. (2003b, JBJ hereafter) obtained an alternative expression for the aperture mass skewness in terms 
of the shear 3PCF and applied this to the CTIO cosmic shear survey, finding a signal at about the 2.3-cr level. 

Following a different approach, the 3PCF was considered directly in a number of recent papers. Schneider & 
Lombardi (2003; hereafter Paper I) defined special combinations of the shear 3PCF which we termed the 'natural 
components', because they obey simple transformation laws under coordinate rotations. In particular, we derived the 
behavior of the 3PCF under parity transformations, and showed that all eight components are expected to be non-zero 
for a general triangle configuration. Zaldarriaga & Scoccimarro (2003) and Takada & Jain (2003a) obtained analytic 
approximations and numerical results, using ray-tracing simulations, for the 3PCF. Takada & Jain (2003b) provide 
an extensive study of the expectation for the shear 3PCF in terms of the halo model of the large-scale structure, and 
they verified the accuracy of their analytical results with numerical simulations. Schneider (2003) investigated the 
transformation properties of a general 3PCF of a polar under parity transformations and showed that the expectation 
value of any quantity containing an odd power of B-modes vanishes for parity- invariant shear fields. 

In this paper, we first consider the relation between the shear 3PCF and the bispectrum of the underlying conver- 
gence (or projected density) field. If the shear field is derivable from a scalar potential (that is, a pure E-mode field), 
as expected for cosmic shear in the absence of intrinsic galaxy alignments and systematics in the observing process, 
the bispectrum of the convergence fully encodes the third-order information of the random field. 1 In this first part 
of the paper, we therefore generalize the relations between the power spectrum of the convergence and the various 
second-order shear statistics (derived in Crittenden et al. 2002; Schneider et al. 2002) to third-order shear statistics. 
After some preliminaries in Sect. 2, we derive the shear 3PCF in terms of the bispectrum of the convergence. From 
these explicit relations, general transformation laws of the 3PCF can be directly seen; for example, the behavior of 
the 3PCF under parity inversion as studied before in Paper I and in Schneider (2003) can be explicitly verified, as will 
be shown in Sect. 4. In Sect. 5 we invert these relations, i.e., we express the bispectrum in terms of the 3PCF of the 
shear. We obtain two formally different expressions for the bispectrum which must, however, be identical in the case 
of a pure E-mode shear field. 

In the second part of this paper (Sect. 6), we consider the third-order aperture mass statistics as a particularly 
convenient integral over the shear 3PCF; in fact, this part of the paper will quite likely be most relevant for future 
studies of higher-order cosmic shear statistics. We first express (M^ p (0)) in terms of the bispectrum and then replace 
the bispectrum in terms of the 3PCF. This procedure yields the same result for (M% (9)) in terms of the shear 3PCF 
as derived by JBJ. We then argue that the third-order aperture measures contain only part of the information about 
the bispectrum of the underlying convergence, and generalize the aperture measures to the case of three different scale 
lengths. We show that this generalization allows us to obtain essentially the full information about the bispectrum. 
These generalized aperture measures are then expressed in terms of the shear 3PCF. Sect. summarizes and discusses 
our results. 

It must be stressed that all our results are valid for other random fields which share the properties of that of cosmic 
shear: A homogeneous, isotropic, parity-symmetric random field of a polar. The most obvious example of such a field 
in the cosmological context, apart from cosmic shear, is the polarization field of the Cosmic Microwave Background 
(e.g., Zaldarriaga et al. 1997) 

2. Preliminaries 

In the first part of this paper (through Sect. |5J) we shall consider a shear field 7 which is caused by an underlying 
projected density (or convergence) field k, as is expected for a shear field produced by light propagation in an inhomo- 

1 It should be noted that even for a pure cosmic shear field, B-modes do occur if source galaxies are clustered (Bernardeau 
1998; Schneider et al. 2002); furthermore, they can occur from slight violations of the lowest-order approximations employed 
when considering light propagation through a slightly inhomogeneous Universe - see, e.g., Bernardeau et al. (1997), SvWJK, 
Jain et al. (2000). 
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Fig. 1. Definitions of the geometry of a trian- 
gle. The X; arc the vertices of the triangle, 
the x; the corresponding sides, tpi are the ori- 
entations of the sides relative to the positive 
Xi-direction, the 4>i are the interior angles of 
the triangle. The figure also shows the ortho- 
center H, i.e. the intersection of the heights of 
the triangle. 



geneous Universe (e.g., Blandford ct al. 1991; Miralda-Escude 1991; Kaiser 1992; see also the reviews by Mellier 1999 
and Bartclmann & Schneider 2001). The relation between the shear (expressed throughout this paper as a complex 
number) and the convergence is most simply given in Fourier space, 

^ i (£) = C ^k(£) 1 (1) 

where (3 is the polar angle of £, and £ is the Fourier transform variable of the angular position vector on the sky. 

In Paper I we considered the 3PCF of the shear. Since the shear is a two-component quantity, the 3PCF has 8 
independent components. Since one cannot form a scalar from the product of three shears alone, one needs to project 
the shear with respect to some reference directions. The three points at which the shear is considered form a triangle, 
and one can project the shear along directions attached to such a triangle, i.e., which rotate with the triangle in 
coordinate rotations. We have considered a number of such obvious projections, namely with respect to the directions 
of the vertices towards one of the centers of a triangle. Let (j be the polar angle of the vector connecting the point X, 
with the chosen center, then the Cartesian components 71.2 of the shear 7 = 71 + i72 are used to define the tangential 
and cross components of the shear relative to the chosen direction Q, 

7(X j; C,) = 7 t(X i ;Ci) + i7x(X i ;C i ) - - [71 (X,) + i 72 (X i )] e" 2i ^ = - 7 (X,)e- 2i <- . (2) 

If the reference directions arc changed, the tangential and cross components of the shear will change correspondingly. 
In particular, defining the 3PCF of the shear components, they will change if a different center of the triangle is chosen. 
In Paper I we defined four complex 'natural components' of the shear 3PCF which show a simple behavior under such 
transformations; they have been termed (x\, x 2 , x 3 ),H = 0,1,2,3. The r«, i = 1,2,3 can be obtained from one 
another by permutations of the arguments Xi , which represent the sides of the triangle. 

We specify our geometry in the same way as in Paper I (see Fig. 1): Let Xj be three points at which the shear is 
considered. The connecting vectors are Xi = X3 — X2, x 2 = Xi — X3, and X3 = X2 — Xi. We denote the polar angle 
of the vector Xj by ipi, and the interior angle of the triangle at the point Xj by fa. Furthermore, we assume that the 
triangle is oriented such that Xi x X2 = X2 x X3 = X3 x xi > 0, where for two two-dimensional vectors a and b we 
defined a x b = a\b 2 — a 2 b\. Hence, the points Xj are ordered counter-clockwise around the triangle. 

In the first part of this paper, we shall consider the projection of the shear relative to the orthoccnter of the 
triangle. As the vector connecting the point Xj and the orthocenter is perpendicular to the vector Xj (see Fig. 1), one 
has d = ipi + 7r/2 in this case, so that 

7 <°>(X i ) = 7t (o) (X0 + i 7 < o) (Xj) = [71 (Xj) + i 72 (Xj)] e" 2i ^ = 7 (X.) e" 2 ^ . (3) 

The shear 3PCF depends linearly on the 3PCF of the convergence, or equivalently on its Fourier transform, the 
bispectrum. The bispectrum of the surface mass density is defined as (see, e.g., van Waerbeke et al. 1999) 

(K(*i)/sCe 2 )«(40) = (2n) 2 [B(£ 1 ,£ 2 ) + B(£ 2 ,£ 3 ) + B(£ 3 ,£ 1 )} §(£ 1+ £ 2 +£ 3 ) , (4) 

hence, it is non-zero only for closed triangles in £-space. This follows from the assumed statistical homogeneity of the 
random field k. Furthermore, if k is an isotropic random field, the function B(£,£') depends only on \£\, \£'\, and the 
angle ip enclosed by £ and £' . We shall therefore write B(£,£') = b(\£\, \£'\, <p). If, in addition, the statistical properties 
of the field k are invariant under parity transformation (as we shall assume throughout), then b is an even function of 
tp, or, equivalently, b is invariant against exchanging £ and £', 



b{i,e,-<p) = b{i,e,<p) = b(t!,i,<p) 



(5) 
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3. The shear three-point correlation function in terms of the bispectrum 

In this section, we will derive the shear 3PCF in terms of the bispectrum B of the convergence. As it turns out, the 
calculations are fairly cumbersome, owing to the rich mathematical structure of the 3PCF with its three arguments, 
compared to the 2PCF which has only one argument (and for which the direction along which the shear components 
are measured is uniquely given by the connecting vector between any pair of points). 



3.1. The case ofT® 

The natural component of the shear 3PCF measured relative to the orthocenter reads 

d 2 h 



rW( Xl ,x 2 ,x 3 ) := ( 7 (o) (X 1 )7 (o) (X 2 )7 (o) (X 3 )) 



d 2 e 2 



(2tt) 2 J (2tt) 2 J (2tt) 



d 2 £ 3 



exphi(^ 1 .X 1 +£ 2 -X 2 +€ 3 .X 3 )] 



exp [2i($>-£ 



9> 



{h(£x)k{£ 2 )h{£ 3 )) 



(6) 



where we made use of the relation (fj^l between the Fourier transforms of the shear and the convergence, and the 
definition ® of the shear components relative to the orthocenter. Inserting the bispectrum 10} into Eq. J§}, performing 
for each of the resulting three terms the integration over the ^-vector which is not in the argument of the B function, 
by making use of the delta-'function' in and using the relations between the corner points and the side vectors 
Xj , one obtains after renaming the dummy integration variables 



T^(x u x 2 ,x 3 ) = 



dli h 



d£ 



2 t-2 



(27T) 



2tt 



d/3i 



2tt 



e i(£ 2 -Xi-£i-X 2 ) , e i(^ 2 -X 3 -^i-X 3 



d0 2 b(£ 1 ,£ 2 ,0 1 -0 2 ) exp 

J^a-Xs-^i-Xi)' 



2i(EA-E 



•■Pt 



(7) 



The angle 3 occurring in Q is the polar angle of the vector £ 3 = —t\ — £ 2 , so that 



cos2/3 3 = 



£l cos 20 x + 1\ cos 20 2 + U x i 2 cos(/3i + 2 ) 



\ti+£i\ A 

We next rename the angles in the following way: 
1 =9 + <p/2; 2 = 9-ip/2- 3 = 9 + , 



sin 20 3 = 



£1 sin 201 + I 2 sin 2[3 2 + 2lil 2 sin(/3i + 2 ) 
\£i+l2\ 2 



(8) 



so that tp is the angle between i\ and £ 2 , as previously defined, and is the angle between the direction of £ 3 and the 
mean of the directions of £\ and e 2 . Since B is independent of 9, a further integration can be carried out in (|7J), using 
/ d0i J d0 2 — J d6 J dip. Owing to the symmetric form of the three terms occurring, only one of the three terms has 
to be calculated explicitly; we shall consider the first term in the following. From the foregoing equations one finds 
that 



cos 20 



{el + e 2 2 )costp + 2e 1 e 2 



sin 20 = 



{el-e 2 2 )sm<p 



\ix+ti\ \tx+t2\" 
Next, we consider the argument of the exponential, 



with \£i +£ 2 \ = el + e 2 2 + 2eie 2 cos ip 



(9) 



£ 2 ■ xi — £\ - x 2 = e 2 x\ cos (9 - ip/2 - ip x ) - e x x 2 cos (9 + (p/2 - ip 2 ) 

= -£ 2 x x sin [9 1 -(<p + <f> 3 ) /2] - t x x 2 sin [9 1 + {ip + <f> 3 ) /2] 



(10) 



where we have defined 9 = 9' + (ipi + ip 2 )/2 and used the fact that tp 2 — <p\ — n — <j> 3 . Therefore, we can write 

t 2 • xi - £ x • x 2 = -As sin(9' + a 3 ) , (11) 

from which one finds, after expanding the trigonometric functions in i|l(J|) and (|llf) . 



A 3 cosa 3 = (£ x x 2 + e 2 x\) cos 



V + 4>3 



A 3 sin a 3 = [t x x 2 — £2^1) sin 



f + 4>3 



a 3 = ^/{e lX2 ) 2 + {e 2Xl ) 2 + 2e x e 2 x x x 2 cos(<^ + ^ 3 ) . 

Finally, we consider the sums over the angles that occur in Q), 



(12) 



^0 i -^<p i =39 + 0-^<p i = 39'+0+(Lp 1 + <p 2 )/2-<p 3 =39'+0 + (<f>i-<fi 2 )/2 



(13) 
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We can now perform the ^-integration of the first term in Q) as follows: 



J J d9 exp [2i(X>~E 



3 i(^ 2 -Xi-£i-X 2 ) _ g 2i/3 e i(0i-0 2 ) / ^/ e 6i9' e -iA 3 sin(e'+Q 3 ) 



2tt 



_ e 2i^ e i(0!-0 2 ) 



/•2:r 

/ dtf e 6i ^- Q3 ) e- iA3 8111,5 = 27re 2i ^e iWl -* 2) e- 6iQ3 J 6 (A 3 ) , 
Jo 



(14) 



where J„ is the Bessel function of the first kind. Therefore, Q becomes 



T^{x l3 X2,x 3 ) = (2tt) 



d^ 2 ^ 2 



2 71 



dtp b(t lt h,<p) e 2 ^ 



o (M ./„ 

where the A4 and a, are obtained from ((12(1 by cyclic permutations of the Xt, x 2 , X3. 



3.2. The case ofT® 

Next, we calculate the natural component 

rM( Xl ,x 2 ,x 3 ) := ( 7 (o) *(X 1 )7 (o) (X 2 )7 (o) (X 3 ))- , r} ^ 2 , 



(15) 



(2tt) 2 y (2tt) 2 

x e 2i(/3 2 +/3 3 -/3i+ V i- V2 - V 3) (k(-£ 1 )k(£ 2 )k(£ 3 )) 



—4 exp [i (*! • Xi - £a • X a - £3 • X 3 ) 

(2tt) 2 



(16) 



where we made use of the fact that k*(£) = k{—£), since k(x) is a real field. Next, we change the integration variable 
£1 —> as a consequence, /3i — ► /?i + 7r, but this does not change the exponential in (|16fl . Inserting the bispectrum 
in the form 10} , and appropriately renaming the dummy integration variables, one finds 



di 2 e 2 



2 k 



d/3i 



2tt 



o (2tt) 2 Jo (2tt)- Jo Jo 

gi^-Xi-^i-XaJ^ifft+fe-^i) , e i(4-X 2 -£i-X 3 ) e 2i(/3i+/32-/3 3 ) 1 e i(£ 2 -X 3 -£i-Xi) e 2i(/3 3 +/3 1 -/3 2 ) 



(17) 



The further calculation proceeds in the same way as in the case of 1^°). Specifically, we employ the change of angular 
integration variables given in © , evaluate the three exponentials containing products of the form Xj • £j using (|10|) , 
111(1 , and their analogous expressions obtained by cyclic permutations of the xi , and calculating the angular sums in 
the exponentials. The final result reads 



r (1) (xi,a, 2 ,x 3 ) = 27T 



dl x 



d^ 2 lo 



(2tt)2 J q (2tt) 2 Jo 



2^ 



d^ 6(4, e 2) ^ 



e i(0i-0 2 +20 3 ) e 2i(/3- v -a 3 ) J 2 (,4 3 ) 



+ e i(0 3 -0 2 ) e -2i( / 3+a 1 ) J ^ A j + Jfa-4 1 -2* 3 ) f jmp+<p-a 3 ) 



(18) 



The expressions for the other two natural components and I^ 3 ) of the shear 3PCF, which are defined in analogy 
to (116(1 by placing the complex conjugation of the shear at point X 2 and X 3 , respectively, are obtained from l|18|) by 
applying the transformation laws given in Paper I, i.e., even permutations of the arguments. 

The resulting expressions 1(15(1 and (|18|) are not only relatively complicated, but their numerical evaluation also is 
quite cumbersome. Recalling that the relation between the 2PCF of the shear and the power spectrum P K (£) of the 
convergence involves a convolution integral over a Bessel function, one should perhaps not be too surprised that in the 
case of third-order statistics there are three such oscillating factors in the transformation between the shear 3PCF and 
the bispectrum. In a future work, we will investigate numerical procedures with which the integration can be carried 
out accurately; first attempts, using Gaussian quadrature for the two ^-integrations and an equidistant grid for the 
^-integration already yielded satisfactory results. Hence, despite the apparent complexity, the foregoing equations can 
be applied in practice. 



4. Transformation laws 



In Paper I the behavior of the natural components under a change of the order of the arguments was discussed, using 
simple geometrical arguments. We shall now consider this behavior explicitly, using the expressions 1(15(1 and ((18(1 . 
First, consider T^°\ 
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Taking an even permutation of the arguments of T^ ) just changes the order of the terms in the integral of (|15f) and 
therefore leaves P" 10 ' unchanged. Taking an odd permutation of the arguments means that two of the arguments arc 
interchanged, e.g., X\ and a; 2 . Interchanging x\ and X2 corresponds to an interchange of <pi and 2 . Using the property 
JSJl, one can also interchange £\ and £ 2 . From © one sees that these changes imply that /3 — > — /3. Furthermore, from 
112(1. one sees that these changes lead to A 3 — > A3, and a 3 — > — a 3 . Together this implies that these transformations 
lead to a complex conjugation of the first term in i|15|) . From the expressions for A; and on obtained from i|12|) by cyclic 
permutations of the Xk, one finds that the above interchanges of x\ and X2 leads to Ai — > A±, Ai — * A%, o>2 — > — ai, 
ol\ — > — a 2 . This then implies that the second term in i|15|) becomes the complex conjugate of the third term, and vice 
versa. Taken together, we see that an odd permutation of the arguments changes T^ ' — > pC )*^ as already argued from 
parity considerations in Paper I. 

Note that an odd permutation of the arguments in geometric terms means that the orientation of the triangle is 
reversed. We shall now show that this is equivalent to replacing all (j>i by —0/. The motivation for this observation comes 
from the fact that for a triangle with odd orientation, the relations [see Eq. (1) of Paper I] between the orientations 
LPi of the x.j and the interior angles 4>i formally yield 4>i 6 [— 7r, 0] (modulo 2ir), whereas for a triangle with even 
orientation, 4>i S [0, 7r] (modulo 2ir). If we apply this transformation, <f>i — > —<f>i, we can change the integration variable 
<p — > — tp in l|15|) . noting from (JSJ) that b(£i, £2, <p) is unaffected by this change. These two changes together then imply 
that /3 — > — /?, Ai — > Ai, and a, — > — ai [see © and respectively]. Hence, all three terms of l(Po|) are transformed 
to their complex conjugates, as was claimed above. Note that this transformation behavior directly implies that T^ ) 
is real if two of its arguments are equal. 

Next one can consider the transformation behavior of T^ 1 ). Cyclic permutations of the arguments transform T^ 1 ) 
into r^ 2 ) and T^ 3 \ yielding the transformation behavior derived in Paper I. Interchanging X2 and x 3 (and thus 4>2 and 
4>z) should yield the complex conjugate of T^ 1 '. Again, we interchange £\ and £2, which then yields (3 — ► —(3, A\ — > A\, 
ai — > ~cni, and so the second term in l |18| l is complex conjugated. Furthermore, these transformations yield A2 — * A3, 
0.2 — > —oi^, A3 —> A2, «3 — » — «2, which shows that the first term in (f 1 81) becomes the complex conjugate of the third, 
and vice versa, so that r*- 1 -* — > r( 1 )* j as was to be shown. This transformation implies that is real if the last two 
arguments are equal. 

Equivalently, we can also let <f>i — > —</>;, and change the integration variable tp — > — tp. This implies A; — > A;, 
a/ — > — cxi, P ~^ ~Pi an d each term in l|18fl is transformed into its complex conjugate. 

5. Bispectrum in terms of the 3PCF 

Recall the situation for the two-point correlation of the shear: there, the relation between the correlation functions 
and the power spectrum of the projected density fluctuations can be inverted (e.g., Schneider et al. 2002; hereafter 
SvWM), and thus the power spectrum can be expressed in terms of the two-point correlation function. It will be shown 
here that in analogy, the bispectrum B{ti,t2) can be expressed in terms of the 3PCF. 

5.1. Bispectrum in terms ofT^ ' 
From one finds that 

(k(£i)k(e 2 )k(£ 3 )) - e" 2i 2> (7(^)7(^)7(4)) 

= e - 2i Eft J d 2 XiJ d 2 X 2 J d 2 X3e 2i E^e i ^- Xl+£2 - X2 +- e3 - X3 )( 7 ( )(X 1 )7( )(X 2 ) 7 ( )(X 3 )) , (19) 

where in the second step we used (01 • We now split the right-hand side of the foregoing equation into three identical 
terms, each of which is thus one third of the above expression, and substitute Xi = X3 + x 2 , X 2 = X 3 — Xi in the 
first of these, and similar substitutions, obtained by cyclic permutations of the X; and x; in the other two terms. This 
then yields 

(k(£i)k(£ 2 )K(e 3 )) = ^- 2i Eft J d 2 Xl J d 2 x 2 J d 2 X 3 e 2i E^e i ( £l - X2 - £2 - Xl )e iX -^+^+ £3 )r( )(x 1 ,x 2 ) + 2 terms , (20) 

where we have used the notation r(°)(xi,x 2 ) for the 3PCF (xi, X2, X3), with x% = |xi + x 2 |, if xi x x 2 > 0, 
and r^>(x 1 ,x 2 ,x 3 ) if xi x X2 < 0. Hence, if expressed in terms of the arguments (xi,x 2 ), the information about 
the orientation of the three points X/ is included. Another notation to be used later on is (xi, X2, ^3), which 
also includes the orientation of the three points. The X 3 -integration in the previous equation yields a delta-function; 
comparison with 10} then results in 

B(£i,£ 2 ) = ie" 2i I> J d 2 xil d 2 x 2 e 2i E^e i A X2 -^ Xl )r( )(x 1 ,x 2 ) . (21) 



P. Schneider, M. Kilbinger & M. Lombardi: The three-point correlation function of cosmic shear II 
It is easy to show that (|21[1 is compatible with 0, since 

r(°)(x 1)X2 ) =|^|^B(£ 1 ^ 2 )exp[2i(^/3 i -^^)] [ e K^-x 1 -£ I .x 3 ) +ei (£,x 2 -£,X3) +e i(£,X3-£,x 1 



(22) 



(£ 2 'Xi-£ r x 2 ) +c i(4-x 2 -£rx 3 ) n i(£ 2 -x 3 -^i-xi 



+ e" 



r (0) (yi,y 2 ) 



where the ip[ are the polar angles of the y i; and y 3 = — y 1 — y 2 . The ^-integrations can be carried out, yielding delta- 
'functions': for the first term, e.g., they yield (27r) 2 <5(y 2 — x 2 )(27r) 2 <5(x 1 — y^), so that also Y^'Pi ~ S Pi- Together, 
this yields 



r<°) ( Xl , x 2 ) = 1 [r(°) (xi, x 2 ) + r(°) (x 2 , x 3 ) + r(°> (x 3 , Xl ) 



(23) 



but since rW(xi,x 2 ) = r(°)(x 2 ,x 3 ) = r(°)(x 3 ,xi), the compatibility of l|21|) and has been shown. 

Next, we want to calculate b(£i,£ 2 ,ip) from (|21|l . i.e. taking a further integration. For that purpose, we write 
fi = /i + tp/2, f 2 = H - VV 2 > where ip x - ip 2 = ip — 3 - Then, 



£i • x 2 - £ 2 • xi = ^ia; 2 cos(# + ip/2 - n + ip/2) - £ 2 x x cos(6> - ip/2 — n — i/j/2) 

= -l\Xi sin \p! - (ip + 3 )/2] - £ 2 Xi sin [// + (<p + 3 )/2] = -A 3 sin(// - a 3 ) , 



(24) 



where we defined // = fi — 6 in the second step, and A 3 and a 3 are given in (|12() . Writing the polar angle of x 3 as 
tps = fx + Cp, one finds that <Pi ~ Pi = 3/x' + ^ — /3, where 



-(a? + £§) cos0 3 + 2a;ia; 2 . (x 2 - x\ ) sm 3 

cos2<^ = ■ -J ; sm2y> = ■ , — 

|xi+x 2 | |xi+x 2 | 

and (3 is given in @. Taken together, l|21|l becomes 

1 



(25) 



dxi x\ / dx 2 x 2 

n Jo JO 



2tt p2ir 



2tt 



~ 2i n dx lXl dx 2 x 2 d0 3 r(°)(x 1 ,x 2 ,0 3 )e 2 ^e 6ia3 J 6 (A 3 ) . 
Jo Jo Jo 



(26) 



It is easy to see that b(£i,£ 2l (p) as given by is real: Since b*(£i,£ 2l ip) = b*(£i,£ 2 ,—ip) - see J5J - and 

r(°)*(xxi x 2 , 3 ) = r(°)(a;i, x 2 , — 3 ), taking the complex conjugate of (ffi))> and simultaneously replacing ip — > — ip, 
and the integration variable 3 — > — 3 yields the same expression as (126(1 . b*(£i,£ 2 ,(p) = b(£i,£ 2 ,<p). 



5.2. Bispectrum in terms of the other 

Next we calculate the bispectrum as a function of the other three natural components of the shear 3PCF, starting 
from 

(Kl(*l)Kl(*2)Kl(4»)> = <Kl(-*l)«l(*2)Kl(40) 



= c -2i(-/5 1 +/3 2 +/3 3 ) 



where we used 



J d 2 X x J d 2 X 2 J ^ X3e i(£ x .X 1 +£ 2 -X 2 +£ 3 -X 3 ) e 2i(- vl +<p 2 + v , 3 ) ^ 7 (o)*(x 1 ) 7 (°)(X 2 ) 7 ^(X 3 )^ (27) 



= J d 2 X e i£ - X 7 (°)*(X) e- 2i¥ 



(28) 



In complete analogy to (|21|l . we split the right-hand side into three terms and make appropriate substitutions. From 
a comparison with (0} we then obtain 



1 d 2 ^! d 2 x 2 e 1 ^'^-^ ^) 



+ I d 2 x 2 y d 2 x 3 S 1 *-**-**-^ + J d 2 Xl J d 2 x 3 e'C^-^-^x-Xs: 



r( 1 )(x 1 ,x 2 )e 2i( - Vl+¥ ' 2+¥,3) . 



(29) 
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Here, we wrote the 3PCF as (Xi)^ (X 2 )j^ {X 3 )) = r( 1 )(x 1 ,x 2 ). When renaming the integration variables, we 
have to apply the transformation rules to the 3PCFs (see Paper I). For the second term, we perform the substitutions 
x 2 — > x 1; x 3 — > x 2 , so that r( 1 ^(x 1 ,x 2 ) — > rW(x3,xi) = r( 3 )(xi,x 2 ). The third term is transformed similarly, and we 
get 

B(£ 1 ,e 2 )+B(£ 2 ,£ 3 )+B(£ 3 ,£ 1 ) = I e -2i(-/3i+/3 2 +&) j ^ j ^ \^i{- V r+^ +93 ) e i{£^-U^) T (i) 



+ e 2i( V i+ V 2-p3) e i(^2-X2-^3-Xi)p(3)/ Xi)X2 -j + e 2i( vl - V2 +V3) c i(^3-X 2 -£i-Xi) r (2)^ x ^ X2 ^ 



(30) 



Unfortunately, the three terms on the right-hand side are not equal, as was the case for i[2()|l. Therefore, we repeat the 
above procedure with (ki(£-l)k*(— £ 2 )kx{£ 3 )) and (ki(^i)ki (€2 )«;*(— ^3)}- The two resulting equations are 

B(£ 1 ,£ 2 )+B(£ 2 ,£ 3 )+B(£ 3 ,£ 1 ) = ~e- 2i ^-f 3 *+^ f d 2 Xl f d 2 x 2 ^(.i-^+.3) e ^rx 2 -l,x l)r (2) (xiiX2) 



+ e 2i(- Vl + ¥ > 2 + V3 ) e i(€2-X 2 -£ 3 -X 1 ) r (l)^ Xi;X2 ^ _j_ e 2i( Vl + V 2-V3) e i(4-X 2 -^i-X 1 ) r (3)/ Xi)X2 \ 



(31) 



and 



B(£i,£ 2 ) + B(^2,4) + S(£3^i) = ie-^^+^-ft' y d 2 ^ d 2 a; 2 p^i+^- V 3) e i(^oc 2 -4ocO r (3) (xij ^ 



g 2i( V i-V2+y3) e i(4-X2-4-Xi)p(2)/ XljX2 ^ + e 2i(- Vl + V2 + V3 ) e i(4-X 2 -£ 1 -X 1 ) r (l)(- x ^ X2 N 



(32) 



Now, we can sum equations 
equal terms, therefore 

Zgypi, P2, p 3 ) 



after moving the /3-phase factors to the left-hand side. Then, we indeed get three 



d 2 zi / d 2 z 2 e 1 



(£i-x 2 -4-Xi) 



e 2i(-vi+P2+V3) r (i)( Xl)X2 ) +e 2i(lpl - V2+¥,3) r^(xi,x 2 ) +c 2i(v)1+tp2 ^ 3) r^(xi,x 2 ) 



(33) 



with 

g{Pl,fa,fa) = e 2i (-^+^+ft) +e 2i(ft-/32+/3s) +e 2i(ft+/32-/3s). (34) 

Equation 133|l can be written as a function of r*- 1 ) only, again using the transformation properties of the 3PCFs 
r( 2 )( Xl ,x 2 ) = r«(x 2 , - Xl - x 2 ) and r( 3 )( Xl ,x 2 ) = r< 1 )(-x 1 - x a ,xi). 

We use H17JI and its counterparts for and r'- 3 -' and insert them into (|33|l . The <^-phase factors cancel, so that 
one obtains 

3 i(k 2 -Xi-k 1 -X 2 ) e 2iA_ ++ + e i(k 2 -X 2 -k 1 -X 3 ) e 2iA + + _ , e i(ka-X s -ki-Xi) e 2iA + _ + 



,i(k 2 -x 2 -k 1 -x 3 ) e 2iA_ ++ , e i(k 2 -x 3 -ki-Xi) ( 



2iAj 



+ c 



^(ks-Xs-ki-Xi) p 2iA_ 



+ e 



i(k 2 -Xi-ki-X 2 )„2iA 



+ + - -I- n 1 



(k 2 -x 1 -k 1 -x 2 ) e 2iA + _ + 

(k 3 -Xa-ki-X 3 )„2iX + _ + 



(35) 



with A±.±.± = ±Ai±A 2 ±A3, where A^ is the polar angle of the vector k^. The x\- and x 2 -integrals can be performed and 
yield 8- functions. This makes the ki- and fc 2 -integrals trivial, yielding k^ = ti and \ = Pi. All the phase exponentials 
add up to give 3g, canceling the pre-factor in (|35|l . leaving only B(£i,£ 2 ) on the right hand side and thus verifying 
C2). 



As was the case for (12611 . one additional angular integral can be carried out. With 



g(Pi,fa,f3 3 ) = c 



210 



(36) 



we find for the three terms: 

g~ 1 e 2i (- , Pi+>P2+'P3) — ^-lg2i(M'+^-03) . ^-lg2i(vi-V2+¥'3) _ g-l e 2i(p,'+ip+4> 3 ) . ^-l e 2i(^>i+(p 2 -(p 3 ) _ --l e 2i(//-^) ^gy-j 

In all three cases, we get the integral J 27r d/x' exp(2LV) exp [— iA 3 sm(p,' — a 3 )] — 27rexp(2ia 3 )J 2 (A 3 ). Finally, 



2 roc />oo />2tt 

Kh,^2,V>) = / dx i x 1 / dx 2 x 2 / < 
3ff Jo Jo Jo 



d</> 3 e^ IQ3 h(A 3 ) 



(38) 
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5.3. Comments 

After having seen the relations of the 3PCF in terms of the bispectrum, one is not surprised to find that their 
inversion derived in this section also is of considerable complexity. This can again be compared to the case of second- 
order statistics, where the power spectrum can be written in terms of the correlation function through an integration. 
The integration extends over all angular scales (as is also the case here) , and so the direct inversion will always be of 
limited accuracy since the correlation functions can only be measured on a finite range of angular scales. In order to 
see the range of application of the previous relations, numerical simulations are probably required. 

The foregoing equations also allow us in principle to express the 3PCF in terms of F^ ', by using the expression 
<|17|) for rW and substituting the bispectrum in this equation by r(°), using 1|21[1 . Whereas the corresponding equations 
can be reduced to a three-dimensional integral, they are fairly complicated; therefore, we shall not reproduce them 
here. 

Given the measured correlation functions from a cosmic shear survey, there is no guarantee that the bispectrum 
estimates from (1261) and l|38[l will agree, even if we ignore noise and measurement errors. The two results will agree only 
if the shear field is derivable from an underlying convergence field, i.e., if the shear is a pure E-mode field. Significant 
differences between the two estimates would then signify that there is a B-mode contribution to the shear. In the 
case of second-order statistics, the separation between E- and B-modes is most conveniently done in terms of the 
aperture statistics (see Crittenden et al. 2002, hereafter CNPT); we shall therefore turn to the aperture measures of 
the third-order shear statistics in the next section. 

At first sight, it may appear surprising that the expression l|26(l for the bispectrum is always real, even though 
we have not constrained to correspond to a pure E-mode field (in fact, we would not really be able to put this 
constraint on the 3PCF - compare the 2PCF: only by combining the two correlation functions £± can one separate 
E- from B-modes). The only assumption we made was that the shear field is parity invariant. This can be understood 
as follows: We can describe a general shear field by the Fourier transform relation Q if we formally replace the 
convergence by n(K) = k e (X) + iK B (X), where k e gives rise to a pure E-mode shear field, and k b corresponds to a 
pure B-mode shear (SvWM). Considering the triple correlator of this complex k, one finds that its real part consists 
of terms ((k e ) 3 ) and (k e (k b ) 2 ), whereas the imaginary part has contributions ((k e ) 2 k b ) and ((k b ) 3 ). As shown by 
Schneider (2003), the latter two terms are strictly zero for a parity- invariant shear field, so that the fact that (|26() is 
real is fully consistent with the vanishing of the imaginary part of the triple correlator of the complex k - both are 
due to the assumed parity invariance. This argument then also implies that the resulting expression Ij26(l for b contains 
both E- and B-modes. One can separate E- and B-modes of the bispectrum by suitably combining the expressions 
11261) and (|38fl . As we shall discuss in Sect.[7| the E-mode bispectrum is obtained by 

6 E = [6(EqOU + 36(EqEHJl] /4 . (39) 



6. Aperture statistics 

We have seen that it is possible to calculate the 3PCF in terms of the bispectrum, and in principle also to invert 
this relation. However, the resulting integrals are very cumbersome to evaluate numerically, owing to the various 
oscillating factors. It therefore would be useful to find some statistics that can be easily calculated in terms of the 
directly measurable 3PCF, but which can also be easily related to the bispectrum. In their very interesting paper, 
JBJ considered the aperture measures, which have been demonstrated to be very useful in the case of second-order 
statistics. The aperture mass centered on the origin of the coordinate system is defined as 

M ap (9) = J d 2 tff/ e (|i?|)/<#) = J d a 0Q e (|0|)7t(tf) , (40) 

where Ug^) is a filter function of characteristic radius 9, the filter function 
2 '* 



Qe($) = ^J W) - W) (41) 

is related to /7e($), and the second equality in l|4T)| is true as long as Ue is a compensated filter, i.e. J r d-dUe('&) — 0, as 
has been shown by Kaiser et al. (1994) and Schneider (1996). jt is the shear component tangent to the center of the 
aperture, i.e., the origin. Hence, 7t(#) + i7 x ($) = — 7 i9* 2 /|i9| 2 , where here and in the following we use the notation 
that a vector x = (x%, x%) can also be represented by a complex number x — x\ + \xi. Hence, "#* 2 /|t?| 2 is nothing but 
the phase factor e _2lc ^, where (j> is the polar angle of t9. 

The aperture mass as a statistics for cosmic shear was introduced by SvWJK who showed that the dispersion 
(M 2 p (#)) of the aperture mass is given as the integral over the power spectrum of the projected mass density k, 
convolved with a filter function which is the square of the Fourier transform of Ue- SvWJK derived this filter function 
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for a family of functions Ug which have a finite support. These filter functions turn out to be quite narrow, so that 
(M% p ($)) provides very localized information about the power spectrum (see also Bartclmann & Schneider 1999). 
Furthermore, SvWJK calculated the skewness of M ap (9) in the frame of second-order perturbation theory for the 
growth of structure. As it turned out, the resulting equations are quite cumbersome, which is in part related to the 
fact that the Fourier transform of the functions Ug chosen contains a Bessel function. 

CNPT suggested an alternative form of the function Ug. When we write Ug(d) = #~ 2 «($/#), then the filter used 
by CNPT is 

«(*) = ^ (l - y) e~* 2/2 ; u(v) = J d 2 x «(|*|) e"? x = £ e-" 2 / 2 ; Qg (tf) = -?L exp . (42) 

Hence, this filter function does not have finite support; this is, however, only a small disadvantage for employing it 
since it cuts off very quickly for distances larger than a few 6. This disadvantage is more than compensated by the 
convenient analytic properties of this filter. 

A further advantage of using aperture measures is that M ap , as calculated from the rightmost expression in (|40|) . 
is sensitive only to an E-mode shear field (see CNPT and SvWM for a discussion of the E/B-mode decomposition of 
shear fields). Hence, defining the complex number 

M(j9):=M ap (fi)+iM x (0) = J d 2 tf Q e (|0|) [ 7t (0) + i 7x (<?)] = - / d 2 tf Q 9 (|0|) 7 (0) r 2 /|tf| 2 , (43) 

M ap (6) vanishes identically for B-modes, whereas M±(8) yields zero for a pure E-mode field. Thus, the aperture 
measures are ideally suited to separating E- and B-modes of the shear. 

CNPT and SvWM have shown that the dispersions <M 2 p (6»)) and (Mj(0)) can be expressed as an integral over 
the two-point correlation functions of the shear. Since the correlation functions are the best measured statistics on real 
data (as they are insensitive to the gaps and holes in the data field), this property allows an easy calculation of the 
aperture dispersions from the data. JBJ showed that the third-order moments of the aperture measures can likewise 
be expressed by the shear 3PCF, and they derived the corresponding relations explicitly - they are remarkably simple. 
The fact that such explicit results can be obtained is tightly related to the choice of the filter function l|42l) ; for a 
filter function with strictly finite support, the resulting expressions are very messy (indeed, we have derived such an 
expression for the filter function used in SvWJK, but it is so complicated that it will most likely be useless for any 
practical work). 



6.1. An alternative derivation of the third-order aperture mass 

We shall here rederive one of the results from JBJ making use of the results obtained in Sect. El the agreement of the 
resulting expression with that of JBJ provides a convenient check for the correctness of the results in Sect. El In a first 
step, we express (M^ p (0)y in terms of the bispectrum. Using the first definition in (|40|l . we find that 

(M! P (6)) = J d 2 tfr U e (\0i\) J d 2 tf 2 [/ e (|tf 2 |) J d 2 ti 3 Ug(\0 3 \) 

When carrying out the ^-integrations, the Fourier transforms Ug are obtained. Inserting the bispectrum in the form 
Q and integrating out the corresponding delta function, one obtains three identical terms. With Ug{£) — u{6£) one 
finds 

d 2 £i f d 2 t 2 



<M a 3 p (0)) = 3 / — -i J J -^B{£ 1 ,£ 2 )u{e\£ 1 \)u{9\£ 2 \)u{e\£ 1 +£ 2 \) (45) 
d£i£i J d£ 2 £ 2 J dip b(ti,£ 3 ,<p) u{6£i) u(6£ 2 ) u (o^t\ + £\ + 21 x l 2 cos^j . (46) 



(27T) 



We shall discuss this result in the next subsection; here, we want to use l|45|l and obtain an explicit equation for \M ap J 
in terms of the 3PCF of the shear. For better comparison with JBJ, we shall slightly change our notation for the 
3PCF. Up to now we have labeled the sides of the triangle formed by the three points by the vectors Xj as defined 
in Sect. El The corresponding natural components of the 3PCF were then denoted by (x\, x 2 , x^) — r(")(xi,X2). 
We shall now define the three points in the form Xi = X3 + y l7 X2 = X3 + y 2 , and then define 

( 7 ( X )(X 1 ) 7 W(X 2 )7 (X) (X 3 )) =f<°>( yi ,y 2 ) ; ( 7 « ( Xl ) 7 «(X 2 ) 7 «*(X 3 )) - f< 3 )( yi , y 2 ) , (47) 
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where the 'x' denotes an arbitrary projection of the shear components, i.e. relative to an arbitrary choice of reference 
directions. The relation between the T and the T follows simply from the definitions of the separation vectors between 
the points Xj, which is y 1 = x 2 , y 2 = — xj, so that 

fi 0) (yi,y 2 ) = fi 0) (x 2 ,- Xl ) = rf >( Xl ,x 2 ) , (48) 

and analogously for the other components of the 3PCF. Now, from combining (|21|l with (|45|) and using our new 
notation for the 3PCF, one finds 

(Mi p (6)) = J dV / d 2 y 2 (y l5 y 2 ) J u(6\*i\) «(W Wi + M) e" 2 ^ 6^+^ , (49) 

where the subscript 'cart' denotes the Cartesian components of the 3PCF. Inserting the Fourier transforms of u 
from and using u(0|A'|) e" 2i/3 ' = u{d\£i\)lf /\h\ 2 = (9 2 £* 2 /2) e'^ 02 ' 2 (where we used, as before, the notation 
£ = £i + \£ 2 ), (EHJ becomes 



! P (9)) = / dv/ d 2 y2 f(°2 t ( yi , y2 ) /0/ W(M) 2 
(l^l 2 + |4| 2 + |4+4| 2 ) 



x exp 



? 2 

— + i (yi • t-\ + y 2 



(50) 



The ^-integrations can now be carried out, by noting that the exponential is just a quadratic function of the integration 
variables, and it is multiplied by a polynomial. The straightforward, but tedious calculation has been carried out with 
Mathematica (Wolfram 1999); the result of the integration then depends on the y { . Substituting the y t in favor of the 
vectors 

2yi-y 2 2y 2 -y 1 y x + y 2 

qi = — 3 — ; q 2 = — 3 — ; q 3 = 3 — . (5i) 

which are the vectors connecting the center of mass of the triangle with its three corners, one can put the result in 
the form 

<K P W) - J-wj-w r -t(yi:y 2 ) g6 ex P ^ ^ J ■ (52) 

Employing now the transformation laws of the natural components of the 3PCF as derived in Paper I, one sees that 
the squares of the complex conjugates of the can be used to obtain the 3PCF with the shear projected along the 
direction towards the center of mass of the triangle, i.e., 

r£L(yi,y 2 )q?q?q? = -fW ( yi ,y 2 ) |^| 2 |g 2 | 2 |g 3 | 2 . (53) 

After this projection of the 3PCF, the integrand depends only on the absolute values of the y 1 and the angle ip between 
them. By carrying out one more integration, one finally obtains 

/ U 3 /m\ 1 f dj/i Vi [ dy 2 V2 f dip - (0) |gi| 2 |g 2 | 2 |<? 3 | 2 / |gi| 2 + |g 2 | 2 + \qa\ 2 \ (t . A s 
(M ap (9)) = -J—J—J— rU n (y u y 2 ,i>) exp ^ — j , (54) 

with 

i2 _ A Vi ~4yiy 2 cos ip + y% 2 _ y\ - A Vl y 2 cos ip + Ay 2 2 y\ + 2 Vl y 2 cost/- + y 2 

ml — g ; ml — g- i ml — g > (ooj 

and thus 

|gi| 2 + |9 2 | 2 + |g 3 | 2 = ^ {yl + yl - 2/12/2 cos ?/>) . (56) 

The result (|54|l agrees with equation (44) of JBJ, after the projection of the 3PCF has been accounted for. It should 
be noted that in the case considered here, where the 3PCF was explicitly obtained in terms of the bispectrum of the 
convergence which, owing to parity invariance, is real (cf. the discussion in Sect. l5~3*)l . the expression 146JI is real, and 
hence H54fl also is real in this case. This can be seen explicitly, since the value of f cen at ip is just the complex conjugate 
one of that at —ip. 

We like to point out that this derivation has shown two interesting aspects: first, the natural component of the 
shear 3PCF arises naturally in this context, confirming the hypothesis of Paper I that this combination of components 
of the correlation functions is indeed useful. Second, the derivation shows that the projection of the shear onto the 
centroid is the most convenient projection in this particular application. 
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6.2. Generalization: Third-order aperture statistics with different filter radii 



How important is the third-order aperture statistics for investigating the third-order statistical properties of the cosmic 
shear? In order to discuss this question, we shall first consider the analogous situation for the second-order statistics. 
There, as mentioned before, the aperture mass dispersion is a filtered version of the power spectrum P K (£) of the 
underlying convergence; for the function Ug considered here, one has 



(M 2 ap (9)) = 



dt£ 

W) 



q4 «4 



(57) 



hence, the filter function relating P K {£) and (M ap (0)) is very narrow, and unless the power spectrum exhibits sharp 
features, the function (M ap (0)) contains basically all the information available for second-order shear statistics (not 
quite - see below). The analogous equation to l|57|) for third-order statistics is given in i|46|) . The function u is very 
narrowly peaked at around £6 <~ 1, and there is one factor of u for each of the three sides of a triangle in ^-space. This 
implies that in the integration of l)46|) the bispectrum is probed only in regions of ^-space where l\ ~ £2 ~ \£\ +£2 1 ~ 1/6. 
Thus, (M 3 p (0)) probes the bispectrum essentially only for equilateral triangles in Fourier space. For this reason, the 
function (M ap (0)) cannot carry the full information of the bispectrum; it merely yields part of this information. 

On the other hand, (jlrJ|l immediately suggests how to improve this situation: if we define the aperture mass statistics 
with three different filter radii 0^, we can probe the bispectrum at wavevectors whose lengths are £i ~ l/0j, and by 
covering a wide range of 0^, one can essentially probe the bispectrum over the full £-space. Indeed, 



(M ap (0 1 )M ap (0 2 )M ap ^ 3 )) 



d 2 h 



B(£ u £ 2 ) u(6 1 \£ 1 \)u(6 2 \£ 2 \)u(6 3 \£ 1 +£ 2 



(2tt) 2 J (2tt) 2 

+^92\£ 1 \)^e 3 \£2\)^9 1 \£ 1 +£2\)+^9 3 \£ 1 \)^e 1 \£2\)u(92\£ 1 +£2\ 



(58) 



(2tt) 3 



Ahh \ d£ 2 £ 2 / d<pb(?i,ta,<p) 



u{6x£x) u(6 2 £ 2 ) u [9^£\ + £%+ 2£ x l 2 cos<A +2 terms 



which illustrates what was said above. Thus, this third-order statistics is expected to be as important for the third- 
order shear statistics as is the aperture mass dispersion for second-order shear statistics. The fact that we do not gain 
additional information by considering different filter scales for the second-order M ap statistics follows from the fact 
that 



<M ap (0i)M a p(0 2 )) 



(2tt) k[ ' 4 6 



40 2 2 




(01 + 91) 



(59) 



Whereas the fact that the mixed correlator can be expressed exactly in terms of the dispersion at an average angle de- 
pends on the special filter function considered here, it nevertheless shows that one does not gain additional information 
when considering the covariance of -M ap . 

We shall now calculate the triple correlator of M ap for three different filter radii in terms of the shear 3PCF, 
essentially using the same method as JBJ. For that, we first calculate the third-order statistics of the complex aperture 
measure M, as defined in (|43|l : 

(M(9 1 )M(9 2 )M(9 3 )) = (M 3 )(6 1 ,9 2 ,6 3 ) = - fd 2 Xt J d 2 X 2 J d 2 X 3 Q fll (|X 1 |)Q e2 (|X 2 |)Q e3 (|X 3 |) 
x < 7 (X 1 ) 7 (X 2 ) 7 (X 3 )) e- 2i ^+^+« , 



(60) 



where the fa are the polar angles of the vectors X«. Writing, as before, Xi = X3 + y 1 , X 2 = X 3 + y 2 , replacing the 
phase factors by e~ 2l( ^ = X* 2 /\Xi\ 2 , and inserting the definitions of the Qg, one obtains 



(M 3 )(9 u 9 2l 9 3 ) = 



-1 



d a ifc / d 2 2 / 2 f(° a ) I ,(y 1 ,y 2 ) 



(4^)3^^404 

d 2 YY* 2 (Y* +y* 1 ) 2 (Y* +y*) 2 exp 



|Y + yil 

29 2 



y 2 l 



29 2 



29 2 



(61) 



where we set for ease of notation the dummy variable X3 = Y. The y-integration is again over the exponential of a 
second-order polynomial in the integration variable, times a polynomial, and thus straightforward to integrate, but 
tedious. Employing Mathematica does most of the job, though its output needed to be further simplified. The result 
is 



(M 3 ) (0i,0 2) 03) = 



S f dyi yi f dy 2 y 2 



24 



02 



92 



(27T) 



Q 6 



(62) 
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where 

e 2 = 



(-9 2 + 20, 



20 2 3 )\qi\ 2 



s = 

(20? 



q2/)2fl2 
i 1 V 2 V 3 



Q6 



20 3 2 )l&| 2 



(201 



20\ 



h = 



+ 01 , (q 2 ~ q 3 )q*l ®l ~ 0. 



29 2 



ml 



6e 2 



68 4 



h = 



(g 3 -gi)g 2 * 01 -01 
29 2 \q 2 \ 2 6 6 2 



fa = 



+ 01 , (ft - q 2 )q*3 o\ - el 



29 2 



93 



6 9 2 



(63) 
(64) 
.(65) 



The choice of the various quantities defined above was made such that S, Z, and the fi are dimensionless, and that 
they become very simple if all 9i are equal. Consider this special case next, i.e., let 8\ — 9 2 = 9 3 = 9. Then, 9 = 8, 



S=l, h = h = h = 1, and Z = (I 



9ij 



192 1 



I© I 



/(28 2 ). Thus, we recover the result Ij54(l in this case which 



was shown above to agree with the result from JBJ. The difference between (|54f) and (|62H is that the former has 
been derived in this paper from the bispectrum of the convergence, and therefore is strictly real, whereas l|62|) has 
been calculated directly in terms of the shear 3PCF and thus applies to arbitrary shear fields, containing both E- and 
B-modes. For reference, we explicitly give the combinations of the qi appearing in the fi above, 



(92 - 93)9* 



3i/2(2yi e 1 ^ - y 2 ) _ (& - M . 3yi(y x - 2y 2 e^) 



I9i I 



4 2/i ~ iV2 cos tp + yl ' 



192 



Vi - 4yi2/ 2 cos ip + Ayl ' 



(gi -92)93 = 3 yl 



2iyiy 2 sinip 



193 1 



yf + 2y 1 y 2 cos^ + yl 



One expects that ^M 3 ) is symmetric with respect to any permutation of its arguments. Indeed, one can show explicitly 
that l|62|) is symmetric with respect to interchanging 8\ and 9 2 . Performing this interchange, changing the variables of 



integration as y\ — ► y 2 , J/2 — ► yi, ip — ► —ip, and making use of the fact that rcen(j/i, 2/2, ip) = ^cen(y 2 , J/i, 27r — one 
finds that these transformations lead to fi —* f2, f2 —* fi, fa f3, and Z is unchanged. To show the symmetry with 
respect to even permutations of the arguments, one needs to employ the symmetry of r^°\x\, x 2: X3) = T^(x 2 , £3, X\), 
and then use either X\ or X 2 as the reference point in the derivation. This then leads to a cyclic permutation of the 
qi and the /j, and thus leaves l)62|) invariant. 

In Fig. |21 we show the latter part of the integrand in (|62|l for the case of three equal apertures (9\ = 8 2 = 8 3 ) and 
for different aperture sizes. Its zeros, if any, are lines of constant y 2 /yi, because the function only depends on the ratio 
of y 2 and 

We next consider the combination of aperture measures 
(M(d 1 )M(8 2 )M*(6 3 )) = (M 2 M*)(9 U 9 2 ;6 3 ) = - J d 2 X x J d 2 X 2 J d 2 X 3 Q fll (|X 1 |)Q e2 (|X 2 |)Q e3 (|X 3 |) 



•,(o) , 



x ( 7 (X 1 ) 7 (X 2 )7*(X 8 )) e-^'+^-W , 



(66) 



where the semicolon in the arguments of (M 2 M*) indicates that this expression is symmetric with respect to inter- 
changing the first two arguments, but not the third one, of course. Using the same conventions for labeling the vertices 
Xj of the triangle as before, we obtain 



(M 2 M*)(6 U 9 2 ;6 S ) = 



-1 



d a yi / d 2 y 2 ff2 t ( yi ,y 2 ) 



(4ir) 3 9f9' : 

d 2 yf 2 (f*+y 1 ) 2 (f*+y*) 2 exp 



yil 



■y 2 \ 



29\ 



28 2 



29 2 



(67) 



After performing the Y- integration and a few manipulations to express the y i in terms of the q i; we obtain 

S f d 2 yi f d 2 y 2 - (3 ) 



(M 2 M*) (6»i 



(2tt) 2 J e 2 J e 2 

1 -*2"*2"2 
1 9l 92 93 r*2 j-*2 i2 



24 e 6 



Jl J 2 J3 



t(yi>y2) 

1 92 I g3 I 

9 e 4 



/r/2/353 



1 

27 



n* 2 

e 2 93 



e 4 e 2 



nn 



where we have defined 
(91 



93 



9J91 

e 4 



92)93 el( 



1 93 



3 9 4 



(68) 



(69) 



and the fi are as before. This form of the equation is easily compared with the result obtained by JBJ, by setting 



9, so that 9 = 6, S = 1, /1 = f 2 = f 3 



93 



1, and Z = (\qi\ 



192 1 



93 



' 12. This then 
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0=(i',l' ; l') = (1',4',1O') 

7£e 7£e Xm 




iogh/i/e] 

Fig. 2. Contours of the integration function ^ ^ ^ /j* 2 /| 2 /jf 2 e~ z - see - as a function of j/ a and y 2 for 
fixed ^ (upper row: ip — 7r/4, lower row: -0 = tt/2). The left-most of the three columns represents the case where all 
three aperture radii are equal. The function scales with the aperture radius. Note that the imaginary part vanishes 
here because of symmetry. The two right columns show the real and imaginary part of the integrand for three different 
filter radii. The contour lines are logarithmically spaced with a factor of 5 between successive lines, starting with 
10~ 10 . Dashed lines correspond to negative values. 



reproduces their Eq. (49), except for a different labeling of the qi (we considered the complex conjugate shear at the 
point X3, whereas JBJ did this at Xi). 

We now employ again the relation between the natural components of the shear 3PCF in the Cartesian reference 
frame and those measured relative to the center of mass of the triangle (see Paper I), 

f?2t(yi,y 2 )9f M = -f^(y 1>ya ) |gi| 2 1? 2 | 2 \h\ 2 , (70) 

and make the corresponding replacements in QfiSjl. after which one more angular integration can be carried out, to 
obtain our final result 



7lJ2/3ff3 + rr loiii 121^1202 S3 ' ci4 I" |2Ci2 J1J2 



J_ \qi\ 2 \fc\ 2 \fa\ 2 f *2 f *2 ,2 
24 Q6 Jl J2 J3 

9 e 4 ' 1 '=" aaa 1 27 Vk"il 2 l92| 2 |g 3 | 2 e 2 aa 1 e 4 |g 3 | 2 e 2 ■' 1 "- ! ' (71) 

which generalizes the result of JBJ for unequal aperture radii. The proof that this last expression is symmetric with 
respect to interchanging 9\ and 9 2 is the same as the one given above. See Fig. [3] for an exemplary plot of the latter 
part of the integrand. 

The product of the four q^s can be written as follows, 

qiMf = ^ {2 [yf + yt + y\y\ (2 cos 20 - 5)] - Vl y 2 [{yl + y 2 2 ) cos0 + 9i(y 2 - y 2 ) sin^] } (72) 
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= {V, V, V) = (V, 4', 10') 

IZe Im IZe Im 




iog[?/i/e] 

Fig. 3. Contours of the latter part of the integrand in 1)71(1 . The plotted function is e~ z times the term in square 
brackets. The upper and lower row correspond to fixed values of ip = 7r/4 and 7r/2 respectively. In the left two columns, 
the real and imaginary part of the function is shown for the three aperture radii being equal. The right two panels 
correspond to three different radii. The contours are the same as in Fig. |5J 



From the two complex triple correlators (M 3 ) and (Af 2 M*), we can now calculate the four real third-order aperture 
statistics, in analogy to what was done in JBJ, 

(M a 3 p ) (0 U 9 2 , 6 3 ) = IZe [(M 2 M*) (0 lf 2 ; 9 3 ) + (M 2 M*) 3 ; 2 ) + (M 2 M*) (0 2 ,9 3 ; 1 ) + (M 3 ) {6 X , 2 ,6 3 )] /4 , 
(M 2 p M ± ) (0 U 2] 9 3 ) = Im [(M 2 M*) (9 1 ,9 3 ;d 2 ) + (M 2 M*) (9 2 ,9 3 ; X ) - (M 2 M*) (9 u 6 2 ; 6 3 ) + <M 3 ) (0 1> 2 ,0 3 )] / 4 , 
<M ap M|) (Of, 2 , 9 3 ) = TZe [(M 2 M*) (0 U 2 ;0 3 ) + (M 2 M*) {0 U 9 3 - 9 2 ) - (M 2 M*) {9 2 ,0 3 ; X ) - (M 3 ) {0i,0^0 3 )] /4 , 

(Ml) (9 U 9 2 , 9 3 ) = Im [(M 2 M*) {0 1 ,0 2 ; 9 3 ) + (M 2 M*) (0 x ,0 3 ; 2 ) + (M 2 M*) (0 2 , 9 3 ; 9{) - <M 3 ) (0,, 9 2 , 3 )] /4 , 

(73) 

with the same notational convention as used before, e.g., (Af ap M^) (9i;02,O 3 ) = (M ap (0i)M±(0 2 )M±(0 3 )) , which is 
symmetric in the last two arguments, as indicated by the semicolon. These four expressions have very different physical 
interpretations. A significant non-zero value of (A/ ap ) indicates that the E-modc of the shear field corresponds to a 
convergence field n which has significant skewness. This is the signal one wants to measure in future cosmic shear 
surveys, and this term contains the information about the underlying cosmic density field, and thus about cosmology. 
A significant non-zero value of (M ap Af^) indicates the presence of a B-mode in the shear field which is correlated with 
the E-mode. Although lensing can generate such a term with small amplitude, by higher-order lensing effects (caused 
by source clustering, violation of the Born approximation in studying light propagation in the Universe, or multiple 
light deflections - see SvWJK for a discussion of these latter effects), these are probably too small to be detectable. 
Therefore, a detection of a (Af ap A/^) most likely will indicate the presence of a 'shear' not coming from lensing, but 
from, e.g., intrinsic alignment of the galaxies (see, e.g., Catelan et al. 2000; Heavens et al. 2000; Crittenden et al. 
2001; Croft & Metzler 2001; Jing 2002). A significant non-zero value of (Mj) indicates that the shear field violates 
parity invariance, as a B-mode shear cannot have odd moments if it is parity-symmetric (Schneider 2003). Finally, a 
significant non-zero value of (M 2 p M±) indicates a parity invariance violation which is correlated with the E-mode shear 
field. Neither of these two latter terms can be explained by cosmic effects which are expected to by parity- invariant, 
but either indicates an underestimate of the statistical errors (coming from the intrinsic ellipticity distribution of 
the sources and from cosmic variance), or the presence of instrumental systematics or artifacts from data reduction 
(cf. the analogous situation for second-order statistics, where a non-zero value of (M ap M±) would indicate significant 
systematics). 
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As stated above, measuring the third-order aperture statistics (through measuring the shear 3PCF and then using 
the foregoing relations) yields essentially all the information about the bispectrum, provided the latter has no sharp 
features in ^-space. The analogous statement for the second-order statistics is not really true: if one considers a cosmic 
shear survey consisting of several unrelated fields of size $ each, one can calculate the aperture dispersion from the 
shear 2PCF for scales, say, Ss $/4. However, the cosmic shear field contains information about the power spectrum 
of the convergence from all scales; in particular, due to the fact that the shear 2PCF is obtained from the power 
spectrum through a filter function which tends to constant for 0£ — > 0, it contains information over the integrated 
power on large scales. Hence, the second-order aperture statistics do not recover the full information about the power 
spectrum contained in the shear 2PCF for a survey of a given size. In order to make better use of the shear data, one 
should take into account a shear measure which contains the large-scale power, such as the top-hat shear dispersion 
on an angular scale comparable to the size of the observed fields, say at 9 = $/4, which can also be obtained in terms 
of the shear 2PCF (CNPT, SvWM). An analogous situation does not exist for the third-order statistics. This can be 
understood intuitively in the following way: Consider again the survey geometry mentioned above, and assume that 
to each of the independent fields a constant shear is added, corresponding to very large-scale power and/or power in 
the bispectrum. The aperture measures will be unable to measure this constant shear, whereas the shear 2PCF will 
be sensitive to it, as will be the top-hat shear dispersion. However, since one cannot form a third-order shear statistics 
which contains the shear only, i.e., without reference directions (such as the direction to the centers of triangles), such 
a constant shear is expected to leave no trace on the shear 3PCF. This can be seen geometrically as follows: consider 
a triangle of points in a constant shear field. Rotation of this triangle by 90 degrees changes the sign of all shear 
components, and thus the triple product changes sign, for which reason a constant shear yields no shear 3PCF. This 
can also be seen algebraically from and H18|l: The occurrence of the Bessel functions, which behave like £ 6 and £ 2 , 
respectively, for small I (at fixed xt) removes all large-scale contributions of the bispectrum in the 3PCF. This fact 
suggests that indeed the third-order aperture measures recover essentially all information about the bispectrum which 
is present in the shear field. 

One might argue that the skewness of the convergence field, top-hat weighted in a circular aperture, is sensitive 
to long wavelength modes, and so third-order statistics on small scales knows about large scales. This is true, and 
may sound like a contradiction to what has been said above. Looking at the second-order statistics first, (k 2 ) (6) is 
sensitive to the power spectrum on all scales i Ss 2tv/9, and it can be expressed by the 2PCF £ + on angular scales < 26 
(CNPT, SvWM). This is due to the fact that the 2PCF of k is the same as that of £ + . Indeed, if one expresses (k 2 ) 
in terms of the resulting convolution kernel has infinite support; hence, (ft 2 ) cannot be expressed through £_ over 
a finite range, because £_ is not sensitive to the power spectrum on large scales. Something analogous happens for 
the three-point statistics. Whereas one can express (k 3 ) in terms of the shear 3PCF (this in fact is easily done, e.g. 
by first expressing (k 3 ) in terms of the bispectrum, and then replacing the bispectrum in terms of the shear 3PCF, 
using the relations in Sect.EJ, the integration range is infinite. One cannot calculate (k 3 ) from the shear over a finite 
region - in fact, the mass-sheet degeneracy does prevent this. Only on very large fields, where the mean of k can be 
set to zero, can one in principle measure (ft 3 ), and that means, one needs information from much larger scales than 
the size of the aperture. 

7. Summary and Discussion 

In this paper we have considered the relation between the 3PCF of the cosmic shear and the bispectrum of the 
underlying convergence field. Explicit expressions for the (natural components of the) shear 3PCF in terms of the 
bispectrum have been derived. These expressions are fairly complicated, and their explicit numerical evaluation non- 
trivial. The transformation properties of the 3PCF under parity reversal can be directly studied using these explicit 
relations and confirm those derived in Paper I by geometrical reasoning. We have then inverted these relations, i.e., 
derived the bispectrum in terms of the shear 3PCF. Two different expressions were obtained, corresponding to the two 
types of natural 3PCF components: one the one hand T(°\ and T^' , i = 1, 2,3 on the other hand. If the shear is due 
to an underlying convergence field, these two expressions should yield the same result for the bispectrum; in general, 
however, if a B-mode contribution is present, these two results will differ. Drawing the analogy to the E/B-mode 
decomposition for the aperture measures in Sect. 16.21 we have conjectured a linear combination of the two expressions 
for the bispectrum which yields the E-mode only. The orthogonal linear combination then yields the cross-bispectrum 
of the E-mode with the square of the B-mode shear. The fact that the bispectrum is real, provided the 3PCF obeys 
parity invariance, reaffirms the result of Schneider (2003) that for a parity-symmetric field, all statistics with an odd 
power of B-modes have to vanish. 

We have then turned to the aperture statistics, using the filter function that was suggested by CNPT and also used 
by JBJ. As a first step we have used the previously derived expressions for the bispectrum in terms of the 3PCF to 
rederive one of the results in JBJ. Then, by considering the third-order aperture statistics in terms of the underlying 
bispectrum we have argued that the third-order aperture statistics with a single filter radius probes the bispectrum 
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only along a one-dimensional cut through its three-dimensional range of definition, namely that of equilateral triangles 
in £-space. Generalizing the aperture statistics to three different filter radii, the full range of the bispectrum can be 
probed, and, in analogy to JBJ, we have derived explicit equations for the generalized third-order aperture statistics 
in terms of the directly measurable shear 3PCF. We showed that using different filter radii did not yield additional 
information in the case of second-order statistics. 

The filter function used in the definition of the aperture measures was that suggested by CNPT. Whereas it does 
not strictly have finite support, this disadvantage compared to the filter function defined in SvWJK is outweighed by 
the convenient algebraic properties it has; these enabled the explicit derivation of fairly simple expressions. 

Let us summarize the features of the aperture statistics which render them so useful as a quantity for characterizing 
cosmic shear (and other polar fields): 

— The aperture measures can be directly calculated in terms of the shear correlation functions. It is the latter that 
can be measured best from a real cosmic shear survey, as they are not affected by the geometry of the survey and 
holes and gaps in the data field. The expressions of the aperture measures in terms of the shear correlation function 
are easy to evaluate by simple sums over the bins for which the correlation functions have been measured. 

— The aperture measures provide very localized information about the underlying power spectrum (in the case of 
second-order statistics) and the bispectrum (for third-order statistics) and therefore contain essentially the full 
information about the properties of the underlying convergence field, unless its power in Fourier space has sharp 
features (which is not expected for a cosmological mass distribution, since there is no sharply defined characteristic 
length scale). 

— One can easily calculate the aperture measures in terms of the power spectrum and the bispectrum, and hence their 
expected dependence on the cosmological parameters can be derived and compared to the measurements. Whereas 
the aperture measures are just one particular way to form integral measures of the shear correlation functions 
- a different integral measure was defined by Bernardeau et al. (2003) and applied to a cosmic shear survey in 
Bernardeau et al. (2002) - it is a particularly convenient one owing to its simple relation to the bispectrum. 

— The aperture measures are the easiest way to separate E- and B-modes of the shear field. Essentially all E/B- 
modc decompositions for the second-order shear statistics have been performed using the aperture measures, and 
we expect that they will play the same role for the third-order statistics. Furthermore, since two of the four 
independent combinations l|73|l of the aperture measures are expected to vanish because of parity invariance, they 
provide a very convenient way to detect remaining systematics in the observing, data reduction and analysis 
process. 

— The aperture statistics are also easily obtained from numerical ray-tracing simulations, since they are defined in 
terms of the underlying convergence in the first place. Hence, in these simulations one can work directly in terms 
of the convergence instead of the more complicated (due to the various components) shear field. 

From the derivation of the 3PCF as a function of the bispectrum, it becomes clear that the definition of the natural 
components have eased the algebra considerably, compared to the case in which one would have tried to calculate its 
individual components (like 7ttx)- Furthermore, the derivation of the third-order aperture statistics directly requires 
the combination of the shear 3PCF provided by the natural components. As discussed in Paper I there are various 
ways to define the natural components of the shear 3PCF, corresponding to the different centers of a triangle. The 
derivation of the aperture measures in terms of the 3PCF has yielded the result that the projection with respect to 
the center-of-mass of a triangle is the most convenient definition (at least in this connection) . 
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